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IGNITION KINETICS OF CARBON MONOXIDE - OXYGEN REACTION* 

by Richard S. Brokaw 
Lewis Research Center 

SUMMARY 

An analytic solution for the ignition kinetics of the dry carbon monoxide - oxygen re- 
action is obtained by assuming an initiation reaction followed by the chain-branching 
steps CO + O - COg and C0 2 + 0 2 - COg + 20 together with chain termination by quench- 
ing of the electronically excited carbon dioxide molecule, C0 2 and the three-body recom- 
bination of carbon monoxide and oxygen atoms. It is concluded that shock -tube induction 
times on reportedly dry carbon monoxide - oxygen mixtures cannot be explained by this 
mechanism. The shock-tube data can be explained by the following scheme, if a water 

k i 

vapor content of 20 parts per million or less is assumed: CO + 0 2 — > C0 2 + O, 
k 5 k 6 k 2' 

O + H 2 0 — 20H, OH + CO — C0 2 + H, H + Og — ► OH + O. The rate constant for the 

initiation reaction was found to be kj = 2. 5x10® exp(-48 000/RT) liters per mole-second. 
The rate constant for the reaction between oxygen atoms and water at 1600° K was esti- 

o q 

mated by two different methods to be 8x10 and 1. 5x10 liters per mole-second. Induc- 
tion periods for mixtures with deliberate additions of hydrogen can be fit by adding the 

k 3» 

steps O + H 2 — ► OH + H and OH + H 2 — ► H 2 0 + H. 

INTRODUCTION 

The explosive reaction of "dry" carbon monoxide and oxygen has been studied by a 
number of investigators starting with Hadman, Thompson, and Hinshelwood (ref. 1) in 
1932. At first glance, one might assume that the kinetics leading to ignition should be 
rather simple and easily elucidated. On the contrary, as recently as 1964, Dickens, 

a|c 
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Dove, and Linnett (ref. 2) wrote that "there is still not general agreement about the 
reaction mechanism. " 

These prior studies have been concerned with explosion limits in static systems, and 
in such systems two factors tend to obscure the elementary kinetics. First, surfaces 
play a role, and the observed limits are affected by prior surface history and treatment; 
there is a "memory" effect (ref. 3). Second, the reaction is greatly accelerated by 
water or hydrogen, and it is difficult indeed to be sure that all traces of hydrogenous 
substances have been excluded. 

Recently Sulzmann, Myers, and Bartle (refs. 4 and 5) have studied the induction 
period preceding rapid carbon dioxide formation in shock-heated carbon monoxide - 
oxygen - argon mixtures. The shock tube is ideal for such studies because surface effects 
are avoided; there is insufficient time for diffusion of chain carriers to and from the 
walls. Furthermore, there is at least the hope that the water -vapor concentration can 
be held low enough so that there are insufficient collisions involving water vapor to affect 
the kinetics. 

In this report, an analytic solution is first developed for the ignition kinetics of the 
dry carbon monoxide - oxygen system. It is found that the experimental data of refer- 
ences 4 and 5 cannot be explained by the mechanism postulated. Next, a solution is ob- 
tained for the kinetics in the presence of traces of water vapor and hydrogen. This 
scheme can be fit to the experimental results by assuming a water -vapor mole fraction of 
about 20 parts per million. This value is larger than the reported contamination of 1 part 
per million; possible rationalizations are discussed. Finally, the experimental data 
(ref. 5) with added hydrogen are adequately explained. 

These solutions are developed by using mathematical techniques that have been suc- 
cessfully applied to hydrogen -oxygen ignition (refs. 6 to 8). The principal assumptions 
are (1) the induction period is isothermal, (2) the reactants are negligibly depleted during 
the induction period, and (3) reactions among chain carriers and between chain carriers 
and products are also negligible. 


SYMBOLS 

A pre- exponential constant, moles liter -1 

b, c Portmanteau constants defined in equation (13) 

C concentration, moles liter - '*' 

AE dissociation energy of reaction I, kcal mole -1 
i rate of initiation reaction, moles liter -1 sec -1 
K equilibrium constant, moles liter -1 
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k specific reaction rate constant 

P steric factor 

R gas constant 

AS entropy change of reaction I 
T absolute temperature, °K 

t time, sec 

x mole fraction 

Z collision number 

\ exponential growth constant (eqs. (11) and (36)), sec 
v reaction frequency, sec - '*’ 

t induction time, sec 

cp roots of equations (8) and (35), sec - * 


IGNITION KINETICS OF DRY CARBON MONOXIDE - OXYGEN SYSTEM 

In this section, a solution is developed for a kinetic scheme comprised of the gas- 
phase reactions suggested by Gordon and Knipe (ref. 3) 


CO + O 5=± CO$ 
k -i 2 


(I) 


CO| + 0 2 -> C0 2 + 20 


(H) 


CO£ + M — ► C0 2 + M 


(HI) 


CO + O + M — *■ C0 2 + M 


(IV) 


together with an initiation step 
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CO + o 2 


co 2 + o 


This reaction scheme involves two chain carriers - oxygen atoms and an electronically 
excited carbon dioxide molecule, CO|. The differential equations governing the growth of 
O and C0 2 concentrations during the induction period are as follows: 

= i - (i^ 1 + ^ 4 )[OJ + + 2i> 2 )[C0 2 ] (1) 


d[CO*] 

— — = ^it°] ~ ( v .i + v 2 + ^3 ) [ C °2^ 
at 


( 2 ) 


Here i s kj[C0][0 2 ] (the initiation rate) and also skj[CO], v_^ =k_ 4 , i^skgfOg], 
i/g s k 3 [M], and s k 4 [CO][M], The v are reaction frequencies and have the dimen- 
sions of reciprocal time. Since the induction period is assumed isothermal and depletion 
of the reactants is neglected, both the initiation rate and the reaction frequencies are 
taken to be constant during the induction period. 

The initiation rate can be eliminated from equation (1) by introducing new variables 


C Q - [O] + 


(y_l + y 2 + i/ 3 )i 


v l v 2 ~ v i v * ~ + v o + v?) 


'rs 


4' -1 


(3) 


and 


'co: 


- [CO*] 




V l u 2 ~ v * ~ + + 


(4) 
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(Cq and C^q* differ appreciably from the actual concentrations [O] and [CC> 2 ] only at the 
2 . 
very start of the reaction. ) Thus i is eliminated, Cqq* replaces [C0 2 ], and Cq re- 

2 

places [O]. A particular solution to the transformed equations is 


0 o = Aq exp(cpt) 


C C0 2 “ 


( 5 ) 
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Substituting equation (5) into the differential equations yields the relations 

+ ^-1 + ^2^CO* = ® 

2 


and 


^i A o 


- (v 


v 2 + 


v z + ^) A co2: = 0 


( 6 ) 


(7) 


A nontrivial solution requires that the determinant of the coefficients 
zero. This leads to the quadratic equation 


of A q 


and A co * be 
2 


<p 2 + 


O'! + v_i + v 2 + ^3 + V £<P ~ 


1^2 


- v l v 2 


~ v 4^-l + v 2 + "3>] = 0 


( 8 ) 


One root of this equation is always negative; the other may be positive if 
v l v 2 ^ v l v 2 + v 4^ v -\ + v 2 + Thus, there is an explosion limit defined by 


k _l + k 2^2-l + k 3t M l 


(9) 


A very similar result has been obtained by Gordon and Knipe (ref. 3) by equating the 
rates of formation and destruction of O and COjjj. 

Call the root of equation (8) that may be positive <p^ and the root that is always 
negative (pg- Next, define X^ s cp^ and Xg - ~<P2' S enera ^ solution is then 


where 


c O = A o, 1 + A o, 2 ex P(- x 2^ 

C CO* = A CO*, 1 + A co|, 2 ex P(-^ 2 ^ 



( 10 ) 


( 11 ) 

( 12 ) 
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and 


b s= v-^ + v_i + i>2 + 

= k 1 [CO] + k 4 + k 2 [0 2 ] + k s [M] + k 4 [CO][M] (13) 


C . »jK 2 - + K 2 + »$) 

- kjkjtco^] - k^jtCOJM] - k 4 [CO][M](k_ 1 + k2[0 2 ] + k s [M]) (14) 


The pre- exponential coefficients in equation (10) are obtained as follows: Initially the 
concentrations of COg and O are zero. Thus, from equations (3), (4), and (10), 


A 0,1 + A 0,2 = ^-1 + v 2 + v 3 ^ 


(15) 


TO 


* 1 + A CO* 2 - V 1 7 


(16) 


and from equation (7), 


^l-^O, i - (^_i + v 2 + v z + X 1 )A cq * x 


"l A 0, 2 ~ ( u -l + v 2 + u 3 ~ X 2^ A C0 2 , 2 


(17) 

(18) 


Equations (15) to (18) can now be solved to obtain 


x 2 ^-1 + v 2 + u 2 + X 1 )A 


°» 1 ~ X^T 


(19) 


X 1 ( v -i + v 2 + v 2~ x 2* 


°> 2 " X^7 


( 20 ) 


A X 2 V 

C °2’ 1 X 1 + X 2 c 


( 21 ) 
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and 



Because the pre- exponential coefficients of equations (10) have been derived (eqs. (19) to 
(22)), equation (23) can be integrated to give 



At this point, it is worthwhile to consider the relative magnitudes of X^ and Xg. From 
equations (13) and (14) it can be shown that 1 > (4c/b^) > -1. Thus it follows that in the 
flammable region (c > O), X^/O^ < 0. 1715, and (Xj/Xg)^ < 0. 0294. Consequently, the 
terms involving exp^Xgt) die out quickly, and furthermore these terms have smaller 
pre -exponential factors. Hence, to a good approximation, 



In attempting to analyze the experimental results of Sulzmann, Myers, and Bartle 
(refs. 4 and 5), it will be assumed that the temperature is high enough so that reaction 
(in), the quenching step, can be neglected, and the pressure is low enough so that chain 
termination by reaction (IV) is also negligible. In this event the growth of COg concen- 
tration with time is given by 


7 



[CO,] (.* 

since 0. 85 < Xg/ ^ 

Figure 1 shows experimental data on the growth of 
time for a mixture of 20 percent carbon monoxide and 10 percent oxygen in argon at 
1597° K and 1. 145 atmospheres (taken from fig. 2 of ref. 4). In figure 1(a), a linear 
plot, the carbon dioxide mole fraction grows in an accelerating fashion for perhaps 400 
microseconds. Next, there is a brief period of linear concentration growth followed by a 
point of inflection and a subsequent decrease in reaction rate. Sulzmann, et al. (ref. 4) 
define the induction period by extrapolating the linear portion of the curve through the in- 
flection point back to zero carbon dioxide concentration. Thus, their induction time cor- 

4 

responds to a carbon dioxide mole fraction of about 5x10 . 

In figure 1(b), the same data are replotted on semilogarithmic coordinates, and it is 


- l) (26) 

carbon dioxide mole fraction with 


65x10 


100x10' 


i-4 



£ 
.4 


•2 ^ 


'.I 

700 


Figure 1. - Growth of carbon dioxide concentration with time (ref. 4). Oxygen, 10 percent; carbon monoxide, 20 percent; pressure, 
1. 145 atmospheres; temperature, 1597° K. 
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clear that the growth of carbon dioxide is exponential, or nearly so, for the first 400 mi- 
croseconds. (The subsequent falloff must be a result either of reactions among chain 
carriers or depletion of the reactants. ) From figure 1(b) and equation (26), it is found 
that = 1. 54x1 0^ per second. Is this exponential growth constant consistent with the 
postulated reaction mechanism? 

p 

Inasmuch as 4c/b < 1, the binominal expansion can be applied to equation (11) with- 
out undue violence (error less than 20 percent), with the result 

H 

If one assumes that the dissociation of COg can be neglected, the rate constants and 
kg can be chosen to give a reasonable (although not perfect) fit to the data of Sulzmann 
and co-workers (refs. 4 and 5). 

Detailed balancing, however, suggests that k ^ is by no means negligible. The 
magnitude of k ^ in comparison with k^ and kg can be determined by estimating the 
equilibrium constant K for the reaction 


kj^tCO^Og] 
kj[CO] + k_ x + kg[Og] 


(27) 


CO* 


‘-1 


co + o 


(- 1 ) 


The entropy change associated with such a process at 1 atmosphere and 1000° to 3000° K 
can be estimated with some confidence by considering analogous dissociations 
(NOg * NO + O, AS/R = 17. 0; SOg * SO + O, AS/R = 17. 3; and COg CO + O, 

AS/R = 18. 1, based on the data of ref. 9). Thus, if it is assumed that AS/R = 17, 

K = 7x10^ exp(-AE/RT) moles per liter, where AE is the dissociation energy of COg. 

If AE = 20 kilocalories, K lg97 o = 130; if AE = 40 kilocalories, K lg97 o = 0.21. This is 
the range of stability for excited carbon dioxide molecules suggested by Clyne and Thrush 
(ref. 10). 

Next solve equation (27) for k^: 

( 28 ) 

k l * k g[Og] kg[Og] 

The last term on the right may be neglected because K > 0. 21 » [CO] = 1. 75x10“ mole 
per liter. Since k^ is positive, 
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k 2 [Q 2 ] > 

K 


[CO] 


8. 8x1 0 6 


Finally, express kg in collision theory form. 


kg = PZ exp 


^2 

RT 


Because the dissociation energy of 0 0 is 7 kilocalories less than that of C0 0 , E„ > AE -7 

n * 11 i — 

kilocalories. At 1597 K, Z ~ 3. 5x10 , so that the steric factor P must be greater 

than 200! This is clearly ridiculous. 

It follows that the experimental ignition -delay data cannot be explained by reactions 
(I) to (IV). 


IGNITION KINETICS OF CARBON MONOXIDE - OXYGEN SYSTEM IN PRESENCE 

OF WATER OR HYDROGEN 

Since Sulzmann, Myers, and Bartle (ref. 4) report that their nominally dry carbon 
monoxide - oxygen mixtures contained about 1 part per million of water vapor, the ki- 
netics in the presence of traces of water must be investigated next. Furthermore, exper- 
iments were performed with small additions of hydrogen (ref. 5), so that the effect of 
hydrogen can profitably be examined as well. 

Once again, assume the initiation step 


i 

CO + Og ^ COg + o 


The product oxygen atoms react with either water or hydrogen 


k 5 

O + HgO 20H (V) 


0 + Hg 



OH + H 


(HI’) 


to generate hydroxyl radicals that oxidize both carbon monoxide and hydrogen 
10 



(VI) 


OH + CO 


COg + H 


and 


OH + H 2 


1 ’ 


h 2 o + h 


(!’) 


Finally, hydroxyl radicals are regenerated through the step 

k 2 t 

H + Og — *• OH + O (H’) 

Reactions (10, (HO, and (HP) are familiar links from the hydrogen - oxygen branched- 
chain scheme. 

The differential equations governing the growth of O, OH, and H concentrations are 
as follows: 



“ ■ (*V + V 5^ C 0 + V 2 ,C H 


(29) 


d^OH 

~ — = ^ 3 * + 2v 5^ C 0 ~ ( v l* + v 6^ C OH + V 2 ,C H 
at 


(30) 


and 


dCj* 

dt 


" ^3 ,C 0 + ( V 1 • + "e^OH 


v 2 ,C H 


(31) 


Here once again a shorthand notation has been introduced so that a k^,[H 2 ], 
v 2' ~ ^2 , t°2i» v 2' ~ ^ 3 '[H 2 ]» ^5 = kgfHgO], and = kg[CO]. Also, the initiation rate 
has been eliminated from the differential equations by introducing new concentration vari- 
ables, 


C Q - [O] 


(32) 


C OH s [OH] + 


(33) 


v V +v 6 
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and 


C H ■ [H] + — (34) 

v 2' 

(Again Cq H and C H differ from [OH] and [H] only early in the indication period. ) As be- 
fore, a particular solution is of the form = Aj exp(<pt), which leads to a cubic equation 
for <p : 

9> 3 + (t' 1 ,+ ^2* + V 2' +V 5 + V Q^ ( P 2 + \^2 ,U 5 + ^ V S ,+ v ^ v V +v q)\ ( P " 2v 2'( v V + v 5^ v V + = 0 


(35) 


Equation (35) has two small roots - one positive and one negative - and one large negative 
root. In searching for the small roots, a little numerical experimentation reveals that 
the first term in equation (35) can be neglected because these roots are small, and the 
third term can be neglected because the amounts of water and hydrogen are small so that 
and are also small. Consequently, 


<P 1 = ~<P 2 - 


2 ^2»( y 3' + v ^( v V + 
v V +v 2' + v v +v 5 +v e 


1/2 


^ k 2 t ^2]( k 3 , t H 23 + k 5[ H 2 < ^j)( k l , ^ H 2J + 

(k i. + k 3 .)[ H 2 ] + k 5 t H 2°] + k 2’[°2.1 + k 6t C °] 


1 1/2 


B X 


(36) 


(Eq. (36) is in error by at most 5 percent for an added hydrogen concentration of 0. 13 
percent. ) Since the coefficient of cp in equation (35) is the negative of the sum of the 
three roots, the large negative root is 


<?3 = -(*!,+ v 2' + V v + v 5 + u 6 ) 


(37) 


The coefficients Aqjj j, Aqjj g, and Aqjj g are found after considerable algebra 
(and neglecting v g, and relative to <pj, cp^, and <pg) 


A OH, 1 


a 1 *V + X i 
2 v 2' v V +v 6 


(38) 
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i 


(39) 


A OH,2 


1 v 2' ~ X 

2 v 2 , ^,+ ^ 6 


and 


A OH, 3 


^ 0 


(40) 


so that 


[OH] = 


[c 

v l» + "6 _ 


cosh(Xt) - 1 + sinh(xt) 
"2’ 


) 


The growth of CC> 2 concentration with time is governed by 

d [ c ° 2 l . 


l + 


dt 


v 6 [OH] 


which can be integrated to give 


[C0 2 ] ^6 


•'l* + ^6 


sinh(Xt) cosh(Xt) - 




v v t 

v V + ^6 


(41) 


(42) 


(43) 


In the absence of added hydrogen and at large times, equation (43) simplifies to 


[C0 2 ] « 


ki[C0][0 2 ] 


- + _l_\exp(Xt) 
X v v ) 


(44) 


which explains the exponential growth of concentration with time shown in figure 1(b). 
Once again, X is determined by the slope of the line. Furthermore, the initiation rate 
can be determined by extrapolating to an intercept, since t> 2 , = k 2f [0 2 ] is known from 
studies of the hydrogen-oxygen reaction; in any event, v^} is small compared with X" . 

Equation (43) can be used to calculate induction times by assuming that the end of the 
induction period can be characterized by some critical concentration of carbon dioxide. 
Figure 1(a) (p. 8) suggests a mole fraction of 5xl0 -4 . Thus, induction times r can be 
obtained by solution of the transcendental equation 
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Induction time, t, usee 



(a) 10 percent oxygen; 20 percent carbon monoxide. (b) 20 percent oxygen; 10 percent carbon monoxide. 



Reciprocaltemperature, 1/T, °K 


(c) 10 percent oxygen; 10 percent carbon monoxide. (d) 5 percent oxygen; 10 percent carbon monoxide. 

Figure 2. - Comparison of computed and experimental induction times as function of reciprxal temperature. (Induction times have been normalized 
to 1 atm, ref. 5. Curves labeled 0 percent H 2 have been calculated by assuming 20 ppm water vapor. ) 
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(i^i f + C0 9 , crit -pm -v r v -t t 

— 1 — = sinh(Xr) + [cosh(XT) - l] + (Xt) (45) 

y 6 k i x C0 x 0 2 P v V V Q 

From equation (36), it follows that X is proportional to the pressure. Hence, the left 
side of equation (45) is pressure independent as are the coefficients x/ v 2 » and v^J U q. 
Consequently, for a given temperature and composition, equation (45) can be solved to 
obtain a value of Xt that is independent of pressure. Hence, the induction times should 
be inversely proportional to the pressure, as is observed experimentally (refs. 4 and 5). 

At high temperatures, the induction periods are short so that sinh(Xt) — Xt, 
cosh(Xt) 1, and equation (43) reduces to 

[C0 2 ] = k.[C0][0 2 ]t (46) 


Under these conditions the induction period is determined largely by the rate of the initia- 
tion reaction, and the growth of carbon dioxide concentration with time should be more 
nearly linear. Such behavior was, in fact, observed (K.G. P. Sulzmann, private commun- 
ication). 

In figure 2, the experimental induction periods of Myers, Sulzmann, and Bartle 
(ref. 5) are compared with values calculated from equation (45) by assuming that 

_•+ = 5x10"^. The rate constant for the initiation reaction has been taken as 
COgj crit 


kj = 2.5xl0 9 exp(-48 000/RT) liters mole -1 sec - * 


(47) 


This rate constant was chosen to fit the delays at the highest temperatures and also the 
linear portion of the curve of concentration growth shown in figure 1(b). This rate con- 
stant is in close agreement with the value 

k^ = (3. 5±1. 6)xl0 9 exp(-51 000±7000/RT) liters mole"* sec" 1 

which Sulzmann, et al. (ref. 4), obtained from measurements of initial slopes of [COg] - 
time histories in the temperature range 2400° to 3000° K. Rate constants for reactions 
(I T ), (H*), and (VI) were taken as recommended by Fristrom and Westenberg (ref. 11). 

No further assumptions were needed to fit the data with additions of 0. Oil and 0. 13 
percent hydrogen. The agreement between theory and experiment shown in figure 2 is 
gratifying indeed and indicates that both the assumed mechanism and values assigned to 
the rate constants are substantially correct. 
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The data without added hydrogen can be fitted by assuming a mole fraction of water 
vapor of 1. 5x10“ ^ and kg as recommended in reference 11. This value is 150 times 
larger than the reported water -vapor content of about 1 part per million. 

From equation (36), however, it follows that the delay depends on the product of kg 
and the water concentration. The value of kg recommended by Fristrom and Westenberg 
(ref. 11) is based on Kaufmann and Del Greco's measurement (ref. 12) of the reverse re- 
action 20H — HgO + O at 310° K and an activation energy assigned by Kaufmann and 
Del Greco. They found the rate of the reaction 20H — HgO + O was about 1/80 of the col- 
lision frequency. If a steric factor of 1 is assumed (so that the factor 1/ 80 is ascribed 
entirely to activation energy) and if the collision frequency is assumed to be proportional 
to T* , one can estimate that kg at 1600° K might be about twentyfold larger than the 
value calculated from Fristrom and Westenberg’s recommendation. This value implies a 
water concentration of 7 parts per million. 

Or, to take another tack, kg can.be estimated by assuming that the rate constant is 
simply the collision frequency times exp(-Eg/RT) and take Eg as the reaction energy 
(15.9 kcal at 1600° K, ref. 9). hi this manner, a rate constant is estimated some 15 
times larger than the Fristrom and Westenberg value. 

There is a third, and perhaps even better way of estimating kg. In figure 1(b), the 
carbon monoxide mole fraction computed from equation (43) is plotted as a function of 
time. Up to about 400 microseconds, theory and experiment are in almost perfect ac- 
cord, but beyond that time the prediction is too high. A reasonable explanation is that 
the water -vapor concentration is depleted - that most of the water has been converted to 
H and OH. Thus, the slopes of a curve faired through the experimental data in figure 1(b) 
can be used as an indication of the square root of the water concentration (see eq. (36)). 
Such estimates, shown in the inset in figure 1(b), must be high, since carbon dioxide is 
produced at a steady rate by reactions (VI) and (II') even after all the water vapor has 
been converted to H and OH. 

The depletion of water vapor during the induction period can be related to the growth 
of carbon dioxide concentration as follows: The rate of destruction of water is 


d[H 2 0] 

dt 


= kg[H 2 0][0] = i — [H 2 Ojsinh(xt) 
X 


(48) 


This expression can be integrated, at least over the region where X is constant, to give 


[h 2 °] c 

[ H 2°] 


k,- k,kJCO][0 9 ] 

= i [cosh(xt) - 1] s J-? _ Jl exp(xt) 

X 2 2x 2 
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Equations (49) and (44) can be combined to eliminate the time and relate the water and 
carbon dioxide concentrations 


In 


[H 2°I 0 

[h 2 0] 


k 5 [co 2 ] 


2 * [G ° 21 


X + 


( 50 ) 


Computed values of the square root of water concentration are also shown in the inset in 

O 

figure 1(b), when kg is assumed to be 8x10 moles per liter second. The upper curve 
was obtained from equation (50) by using the experimental carbon dioxide concentrations. 
(This curve is clearly an upper bound since COg formation continues by reactions (11*) 
and (VI) even after the water is consumed. ) The lower curve is calculated from the CO 2 
concentration according to equation (44) and is a lower bound. This value of kg implies 
a water -vapor content of about 30 parts per million. 

Finally, after all the water has been converted into H and OH, the carbon dioxide 
concentration should grow in a linear fashion according to 


d[co 2 ] 

dt 


= *i[CO][0 2 ] 


. 2 *2.»WO] 

* 2 ,[° 2 ] + k 6 [CO] 2 O 


(51) 


where [H«0] is the initial moisture content. From figure 1(a), at large times (dashed 
O 

line), d[COg]/dt ~ 0. 16 mole per liter-second, which is large compared with the initiation 

rate. Since values for kg, and kg are available, the initial water concentration can 

now be calculated from equation (48); it turns out to be about 17 parts per million, which 

leads to kp. = 1. 5x10 liters per mole-second. 

0 8 

The estimate of kg = 8x10 (from the dropoff in COg growth rate) can be combined 
with Kaufmann and Del Greco’s value at 310° K to obtain 


kg = 3.8x10^ exp(-19 500/RT) liters mole"* sec"* (52) 

The induction periods without hydrogen addition in figure 2 have been computed by as- 
suming a water vapor concentration of 20 parts per million and kg from equation (52). 
The data can be fit equally well by assuming lower water -vapor concentrations with cor- 
respondingly larger values of kg. The agreement between theory and experiment is good, 
and strongly supports the kinetic scheme assumed provided that one accepts the idea that 
the water -vapor concentration may have been 7 to 30 times the value reported in refer - 
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ence 4. This is perhaps a reflection of the problems that often plague accurate measure- 
ment of small water -vapor concentrations. 


CONCLUDING REMARKS 

The arguments of the preceding sections lead to the following conclusions regarding 
the oxidation of carbon monoxide: 

1. The rate of the reaction CO + 0 2 — COg + O seems well established, with the rate 
constant given by equation (47). The rate of the reverse reaction k_- can also be ob- 
tained by invoking detailed balancing and appropriate thermochemical data (ref. 9): 

k . = 2x10® exp(-53 000/RT) liters mole - * sec - * 

2. The rate of the reaction O + H 9 0 — 20H has been estimated by two different 

O Q " O 

methods as 8x10 and 1. 5x10 liters per mole-second near 1600 K. 

3. The extraordinary sensitivity of the carbon monoxide oxidation to traces of hydro- 
gen and water vapor is revealed once again. Thus, water -vapor concentrations of about 
20 parts per million seem to dominate the kinetics in experiments lasting a few tenths of 
a millisecond. 

4. If parts per million of water vapor dominate the kinetics in the millisecond time 
scale, it seems reasonable that parts per billion of hydrogenous material may well de- 
termine the kinetics in static experiments at lower temperatures where reaction times 
are several seconds or minutes. 

5. The data of Sulzmann, et al. (refs. 4 and 5) cannot be explained by the chain - 
branching-reaction scheme 


co + o 



co| + o 2 



C0 2 + 20 


(I) 

(n) 


This suggests (together with item 4) that there is little or no experimental evidence that 
reaction (n) occurs. 

6. Finally, what can be postulated as to the explosion of absolutely hydrogen-free 
mixtures of carbon monoxide and oxygen? A reasonable guess is that these should be 
thermal explosions, initiated through the nearly thermoneutral reaction 
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co 2 + o 


CO + o 2 


k i 


with subsequent heat release due to 

O + CO + M — C0 2 + M 


and 


O + O + M — Og + M 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, May 23, 1966, 
129-01-02-01-22. 
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